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The optical conductivity ct(cj) and the d.c. resistivity p(T) within the extended t-J model on a square lattice, 
as relevant to high-T c cuprates, are reinvestigated using the exact-diagonalization method for small systems, 
improved by performing a twisted boundary condition averaging. The influence of the next-nearest-neighbor 
hopping t' is also considered. The behaviour of results at intermediate doping is consistent with a marginal- 
Fermi-liquid scenario and in the case of t' = for lu > T follows the power law a <x ui~ v with v ~ 0.65 
consistent with experiments. At low doping < 0.1 for T < J o~(w) develops a shoulder at uj ~ uj* , 
consistent with the observed mid-infrared peak in experiments, accompanied by a shallow dip for u < uj* . 
This region is characterized by the resistivity saturation, whereas a more coherent transport appears at T < T* 
producing a more pronounced decrease in p(T). The behavior of the normalized resistivity cnp(T) is within a 
factor of 2 quantitatively consistent with experiments in cuprates. 

PACS numbers: 71.27.-W, 72.10.-d, 74.72.-h 



I. INTRODUCTION 

Since the discovery of high-T c cuprates, strong correlations 
among electrons have been considered as the crucial reason 
for anomalous transport properties in the normal state of these 
materials. Prominent example is the d.c. resistivity obeying 
the well-known linear law, p oc T, in the intermediate (op- 
timum) range of hole doping accompanied by an anomalous 
but universal frequency dependent optical conductivity a(u>), 
phenomenologically described via a marginal - Fermi - liquid 
(MFL) scenario 1 or the quantum-critical behavior. 2 In the last 
decade the emphasis has been centered on the experimental 
investigations of the underdoped regime in various cuprates. 
The main signatures of this regime are the kink in d.c. p(T) 
at the crossover scale T* 3,4 and the appearance of a broad 
peak in the infrared region (mid-IR) 5 as well as of the pseudo- 
gap scale in <j(oj)&L Recently, it has been established within 
this regime that the mobility of holes increases with doping 
concentration^ 

From the theoretical point of view, it is not an easy task to 
construct an analytical transport theory for strongly correlated 
systems, starting from a microscopic model that describes the 
motion of charge in a weakly or moderately doped antiferro- 
magnetic (AFM) insulator. Even at present such a theory is 
lacking. However, numerical investigations of the prototype 
models, in particular of the t-J model, show that a(iv) and 
p{T) of cuprates can be reasonably accounted for. Applica- 
tion of the finite-temperature Lanczos method (FTLM)2ii£ to 
small systems yields results being most reliable in the regime 
of intermediate dopin g 11 ! 12 where also the experimental win- 
dow T < 1000 K can be probed best. In the latter regime, 
numerical results show an overall agreement with the simple 
scaling er(w) oc (1 — exp(—oj/T))/oj, with the only scale 
given by temperature TpiSi^ 

Recently, several aspects of the conductivity in cuprates 
have been reopened urging for reconsideration and im- 
provement of theoretical and model results. Detailed 
experimental studies in the optimally doped cuprate 



Bi2Sr 2 Cao.92Yo.o8Cu20 8+< 5 (BSCCO)* reveal for u) > T 
the T-independent power-law behavior <t(u>) oc uj~ v with 
v ~ 0.65. The authors attribute such a scaling to the vicin- 
ity of the quantum-critical points 

Another question is the existence and value of the resistivity 
saturation in cuprates at low hole dopingpi^ii It has been re- 
alized that cuprates at high T and low doping show very large 
resistivity p2& and thus violate the naive Ioffe-Regel condition 
for metals, which implies a saturation of p(T) when the scat- 
tering length L s reaches the intercell distance. On the basis 
of the t-J model the modified saturation value p sat has been 
derived 13 and its anomalous large value has been ascribed to 
the kinetic energy being strongly suppressed due to strong 
correlations, in particular (Hkin) oc Ch close to half-filling, 
where c/j = Nh /N is the hole number concentration. From 
the experimental view, La 2 - xStxCuO^ (LSCO) serves as a 
prototype cuprate with a well controlled hole doping = x 
and stable material properties. Previous and recent results^ in 
single-crystal LSCO have revealed signatures of saturation at 
high T > T* ~ 500 K. Plausibly, the saturation phenomenon 
is related to the emergence of a T window, where a(uj < uj*) 
is quite flat or even develops a shallow dip. 1516 The latter phe- 
nomenon is clearly related to the appearance of a broad mid- 
IR peak at u> ~ uj* first found at low doping in LSCC 5 J£ and 
recently also in YBa 2 Cu 3 2/ (YBCO)-i& 

Recently, estimates were presented 1314 that the t-J model 
yields substantially smaller resistivity p(T) than experiments 
in LSCO for the regime of low doping. This apparently leads 
to a conclusion that coupling to additional degrees of free- 
dom, in particular to phonons, might be essential to explain 
the large resistivity. This aspect is clearly of importance since 
it is related to the central question, i.e., to what extent the t-J 
model and strong correlations alone can describe the physics 
of high-T" c cuprates. 

The existence of the mid-IR resonance in <j(uj) is exper- 
imentally well established at low dopingriia*i* however the 
consensus on its origin has not been reached yet. We note 
that such a peak has been reproduced already in calculations 
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of a(uo) for a single hole within the t-J model at T = C 17 
and has been attributed to the string-picture of incoherent hole 
motion leading to the shoulder at lo* ~ 2 J. 18 However, a con- 
firmation of this feature in more realistic T > calculation 
was missing so far. 

The aim of this paper is to give answers to above ques- 
tions via a systematic reinvestigation of (t(lj) and p{T) within 
the extended t-J model, considering also the influence of the 
next-nearest-neighbor (NNN) hopping t' which has already 
been invoked in the modeling of hole and electron doped 
cuprates 19 to possibly reveal the pronounced difference of 
<j(uj) between both classes of cuprates. Additional hopping 
parameter t! has been identified to be important for the ex- 
planation of spectral properties, e.g., the band dispersion as 
measured via the angle-resolved photoemission experiments 
(ARPES). 20 Even more, t' < could be the essential param- 
eter discriminating various families of hole-doped cuprates 
(LSCO, BSCCO etc.), 21 - 22 in particular their superconducting 
transition temperature T c . Relative to the previous workiMi 
we are able to study somewhat larger size systems using the 
FTLM. At the same time we improve the method by introduc- 
ing averaging over twisted boundary conditions (TBC), de- 
scribed furtheron. The improvement shows up in a more con- 
trolled behavior of cr(u>) at low u, being essential to extract 
reproducible p(T) at lower T. In the intermediate doping we 
present calculations for systems up to N = 20 sites to clarify 
the universal behavior of a(oj). In the low-doping regime we 
study systems up to N = 26 sites allowing us to establish the 
emergence of a pseudogap scale in u(lo) as well as the onset 
of a more coherent transport for T < T*. Our results still re- 
main restricted to T > 400 K. Nevertheless, they confirm the 
qualitative behavior of <j{u>) in cuprates, with the quantitative 
discrepancy in p(T) compared to experiments within a factor 
of 2 at most. 

The paper is organized as follows. In Sec. II the improved 
FTLM method, employing the TBC averaging, is described. 
The comparison of results with the usual fixed-boundary con- 
dition (BC) method is presented. In Sec. Ill the optimum dop- 
ing results for cr(uS) and p(T) are presented and discussed in 
connection with the scaling behavior and experimental results 
for cuprates. In Sec. IV the low-doping regime is examined, 
with the emphasis on the emergence of the shoulder corre- 
sponding to the mid-IR peak in a{uS), the resistivity saturation 
and the onset of a coherent transport for T < T*. Conclusions 
are given in Sec. V. 



II. NUMERICAL METHOD 

In the following we study numerically the extended t-J 
model 

H = -y^Jijc\ s Cj s + J ^(S,; • Sj - jriirij), (1) 

«>J> (ij) 

on a square lattice, whereby we include besides the nearest 
neighbor (NN) hopping tij = t also NNN hopping = t'. 
Strong correlations among electrons are incorporated via pro- 
jected operators, e.g., c\ s = (1 — 7i» _«)ct , which do not allow 



for a double occupancy of sites. The dependence of cr(u>) on t' 
has already been examined in connection with the difference 
between the hole-doped and electron-doped cuprates. 19 How- 
ever, the influence of t' on the d.c. transport is less evident 
and has not been studied systematically so farJ Plausibly, 
it is expected that t' < 0, being appropriate for hole-doped 
cuprates, leads to a larger frustration of the AFM spin back- 
ground and consequently to a larger resistivity. In our study 
we will test the behavior at three different t' = 0, t' = —0.15 1 
and t' = — 0.3 t, respectively. We keep everywhere J = 0.3 t 
as appropriate for cuprates, as well as we note that t ~ 0.4 eV 
for a direct comparison with experiments. 

The main limitation to the validity of numerical results 
comes from the finite-size effects which begin to dominate 
results at low T < T/ s . As a criterion for Tf s we use the 
thermodynamic sumlS 



Z(T) = Trexp[-{H-E )/T], 



(2) 



calculated in a given system for fixed hole number Nh and 
the requirement Z(Tf s ) = Z* ^> 1. In the following we use 

Z* ~ 30. 

Finite size effects can be substantially reduced by employ- 
ing the TBC averaging! 2 ^! In a system with periodic BC the 
latter is achieved by introducing an uniform vector potential 8 
modifying the hopping elements iy — > tij — exp(id ■ ry) 
by a phase factor. We use furtheron N t different phases in- 
stead of an usually fixed BC 8 — 0. For an arbitrary operator 
A we then perform the usual thermodynamic averaging over 
N t different TBC 



(A) 



E 



1 Z 3 



(3) 



where (A)j and Zj refer to the canonical FTLM expectation 

values obtained for each fixed phase 8. In order to preserve the 
translational invariance of the lattice Hamiltonian in 2D we 
are allowed to study tilted square lattices with sites where 



N must be Pythagorean, i.e., N 



-m . Consequently the 



first Brillouin zone, where the phases are chosen equidistantly, 
is a tilted square and Nt is also chosen Pythagorean, so that 
the phases form a regular square lattice as well. Effectively, 
such a choice reproduces, e.g., for free fermions the regular 
(square) lattice of NN t points in the k space within the first 
Brillouin zone. 

There are several advantages of the TBC averaging proce- 
dure: 

a) In the limit T — > (for large enough N t ) this method 
reproduces for an arbitrary quantity, e.g. for <r(w), the result 
corresponding to the ground state wavefunction |^o(^)) with 
Eq(6) — min. for a chosen system with N sites and Nh holes. 

b) In a finite (nonintegrable) system, even for T > 0, we 
generally expect the 2D optical conductivity of the form 



where 



a(uj) = 2itD c 5(uj) + a reg (uj), 

1 1 f°° 

2N Treg J 0+ 



(4) 
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Figure 1 : Dynamical 2D mobility a / for a single hole Nh — 1 in 
different systems N = 16, 18, 20, 26 evaluated at fixed T/t = 0.15: 
a) for fixed BC 9 = 0, and b) with the TBC averaging. 



D c is the charge stiffness representing the nondissipative part 
of conductivity, emerging due to nonscattered coherent charge 
propagation in a finite system with periodic BC, and (r) is the 
kinetic stress tensor or the generalized kinetic energy which 
in the case of NN hopping on a square lattice reduces to 
the usual kinetic energy, i.e., (r) = — (Hf-in) /2. 25 Since 



D C {T = 0) oc d 2 Eo(0)/d 2 9\ 



Eq(6)— min. 



> 



the TBC 



averaging leads automatically to D C (T — > 0) > as appro- 
priate for a metal. Note that previous calculations of <j(ui) at 
fixed 9 = 1011 generally suffered on a quite system depen- 
dent D C (T — ► 0) which was in many cases even negative. 

c) In general, the TBC averaging reduces finite-size effects, 
although it cannot completely eliminate them. Results for 
<j(uj) are thus more size and shape independent. In Fig. 1 we 
present the comparison of results for a single hole Nh = 1 
in different systems N — 16, 18, 20, 26, as obtained for a 
fixed BC 9 = and when using the TBC averaging. There 
is a clear difference in the 1ow-cj regime where D c can be 
even negative for a fixed BC, while rather consistent results 
are obtained when applying the improved method. Moreover, 
the TBC averaged u(ui) display quite evident and consistent 
shoulder at lu* ~ 0.5 t, corresponding to the mid-IR peak in 
cupratespii^i^ while this feature is hardly visible for fixed BC 
spectra. 

In the following, we present results using typically N t ~ 10 
and M ~ 140 Lanczos steps within each symmetry sector of 
the Hilbert space. Since TBC averaging now also takes the 
role of random sampling we do not perform any additional 
sampling over the initial Lanczos wavefunction within each 
sector. In each of the available systems N = 18, 20, 26 we 
are able to reach Tf s ~ 0.1 t for the intermediate doping 
and somewhat higher for low doping. In the analysis of a(ui) 



spectra broadening e = 0.07 t is used. 



III. INTERMEDIATE - OPTIMUM DOPING 

One of the most striking facts in cuprates, recognized 
from the beginning, is the universality of the charge response 
a(uj,T) in the vicinity of the optimum doping. Whereas 
the d.c. resistivity shows a linear variation p(T) oc T, also 
<j{lS) can be analyzed with an anomalous Drude-like relax- 
ation rate 1/t(u>,T) oc lo + £T. The only relevant u> scale 
seems to be given by T itself, as summarized within the 
phenomenological MFL scenario. 1 Recent study of optimally 
doped BSCCO system has revealed, that at low T in a broad 
range of T < lo < 0.5 eV one can describe results with a 
power law <j(u>) oc lo~ v , whereby v ~ 0.65. 2 Still, it is not 
possible to represent results for cr(u>, T) in the whole (u>, T) 
range with a single universal function of uj/T as, e.g., required 
within the quantum-critical-point scenario. 2 
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Figure 2: 2D optical conductivity a{uS) for the intermediate doping 
Ch = 3/20 as calculated for N — 20 square lattice for different T/t 
andt'/i = 0, -0.15, -0.3. 



Numerical results for a(ui) within the t-J model are most 
reliable at intermediate doping, since there the lowest Tf s can 
be reached within a fixed system of N sites. 10 In Fig. 2 we 
present results for a lattice of N — 20 sites and Nh — 3 holes, 
i.e., the hole doping c/j = Nh/N = 0.15, at various T/t and 
different t' /t = 0, —0.15, —0.3. Obtained 2D cr(uS) can serve 
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as a test for the scaling behavior. The optical sum rule, 



(6) 



requires a fast fall-off of <j(uj) for large uj > oj c , hence one 
can discuss the scaling only for u> < uj c ~ 2t. 

In Fig. 3 we study a{uj) for t' = within the regime T < 
uj < uj c . We see that the fall-off is slower than oc 1/uj and 
is well described with the power law a oc uj~ v with v < 
1. Still, v seems to depend slightly on t'. The closest fit for 
t' = yields v = 0.65 which is in excellent agreement with 
experiments in BSCCO, 2 although for the latter cuprates more 
appropriate model should be the one with t' < 0. For t' < 
our results reveal a decreasing v. I.e., we get for t' /t = —0.15 
and t' /t = —0.3, v = 0.5 and v = 0.42, respectively. 
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Figure 3: Log-log plot of a vs. uj for = 3/20. Full line represents 
the closest fit to the power law with v = 0.65 for T/t — 0.1. 

In order to come closer to the understanding of the scaling 
behavior, we represent a(uj) in a general form following the 
Kubo formula, 



1 



■j/T 



a{uj) = C{uj) 
where C{oS) is current-current correlation function 

C(w)=Re/ dte^(j(t)j). 



(7) 



(8) 



It has been observed that the scaling behavior of a(uj,T) 
within the t-J model at intermediate doping can be well de- 
scribed up to the cutoff uj < uj c by a parameter-free form of 
Eq. |7} with C(uj) ~ C where C is T-independent. Such 
a form is clearly a restricted version of the MFL behavior 1 
and reproduces both p oc T law and a(uj 3> T) oc 1/uj. In 
view of present improved results presented in Fig. 3 we can 
analyze the deviation from the simple C{uj) — Co form. In 
Fig. 4 we show G(uj) for t' — and various T/t. In this 
case the spectra are broadened with the factor e = O.lt. The 
main characteristic is that C(uj) is an increasing function for 
uj < uj c , consistent with effective v < 1. 

The origin of a very broad range of validity of a oc uj~ v has 
not become clear yet. Note that v < 1 requires an increasing 



C(uj) which in normal metals is supposed to have a Lorentzian 
form with a characteristic width determined by a Drude re- 
laxation rate 1/r < T. The interpretation of the anomalous 
(constant or even increasing) C(uj) has been proposed by one 
of the authors. 27 By performing a simple decoupling of C(uj) 
in terms of the single-electron spectral functions A(k, uj) and 
neglecting the vertex corrections, one gets 



C(uj) 



2ire 



duj 1 x 



f(-uj')f(uj' - uj)A(k, uj')A(k, uj'-uj), 



(9) 



where / is the Fermi function and unrenormalized band 
velocities. We represent A(k, uj) for quasiparticles close to 
the Fermi energy as 



1 



A(k,uj) = - 



k 1 k 



^- £k ) 2 + r2' 



(10) 



where parameters Z^, Tk, £k in general dependent on T. In 
order to reproduce the MFL scaling of <j(uj) we have to as- 
sume the MFL form for the damping, i.e., T — j(\uj\ + £T) 
as well as neglect the k dependence of T and Z. With these 
simplifications we arrive at 



C(w) = C 



ur + 1 (uj,uj') z 



where T(uj,uj') = T(uj') + T(uj' - uj). It has been noted 27 
that C(uj) ~ Co appears already for 7 ~ 0.3, while for 
7 > 0.3 it is an increasing function of uj/T. We should also 
remind that 7 > 0.3 is consistent with ARPES experiments in 
BSCCO, where along the nodal direction 28 as well as in other 
directions 29 7 ~ 0.7 was obtained. 
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Figure 4: C'(uj) for intermediate doping Ch = 3/20 and various T/t. 
Full line represents the analytical result from Eq. lilt with 7 = 1.2, 
£ = yrand T/t = 0.1. 

For comparison with the numerical result at the lowest 
T/t = 0.1, we display in Fig. 4 C(u>) following from Eq. il 11 
with 7 = 1.2 and £ = it. The qualitative behavior is satisfy- 
ing, with a visible deviation only at low uj < T. This indicates 
that anomalous C(uj) as well as the power law <j(uj) oc uj~ v 



with v < 1 can be qualitatively described with the assumption 
that the quasiparticle damping is of the MFL-form with large 
7 > 0.3. Still, the description with the simple Eq. Jl li cannot 
quantitatively account for the whole region (uj, T) < u> c . 
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Figure 5: Normalized resistivity c; l p(T) at intermediate doping for 
t' = and t'/t — —0.15, compared with the experimental data 
for LSCO. Squares denote results for x = 0.15, 3 while circles and 
triangles are for x = 0.15, 0.18, 8 respectively. 

Let us discuss the quantitative comparison of our novel re- 
sults for p(T) = c/a(uj = 0, T) (we use c = 6.6A) with 
the experiments at intermediate doping. In Fig. 5 we com- 
pare experimental results for the normalized resistivity Chp of 
LSCO forx = 0.15, 0.18^ with corresponding FTLM results 
for Nh = 3 in systems with N = 18, 20 sites in the case of 
t' = and t'/t = —0.15, respectively. Results for t' = are 
well consistent with the linear p oc T law in the T regime of 
experimental relevance. It is characteristic that t'/t = —0.15 
results reveal somewhat larger p at lower T while also the 
presumed linearity is reached at lower T < 1000 K. Taking 
into account the variation of different experimental data the 
quantitative agreement is satisfying although we allow for the 
possibility that our model results underestimate p. 



IV. LOW-DOPING REGIME 

In previous numerical studies of cr(uj, T) results for the low- 
doping regime 1011 were less conclusive. The reasons were the 
following: a) Tf s increases as Ch — > 0. So even Tf s ~ 0.15 t 
was unreachable in the systems available with the FTLM a 
decade ago. b) At the lowest doping, i.e., for the case of a 
single hole Nh = 1, there is a significant nondissipative con- 
tribution D c at T ~ Tf s . This indicates that the effective 
scattering length L s might become larger than the system size 
y/~N. In addition, D c is varying quite uncontrollably between 
different system sizes when the FTLM is performed with a 
fixed BC. In fact, as shown in Fig. 1 in several cases one gets 
unphysical D c < 0. c) Systematic experimental studies of 
p(T) and u(lo) in cuprates at low doping have become avail- 
able more recently*^*^ 

The TBC averaging substantially improves the FTLM re- 
sults for a(ui), as seen in Fig. 1. They are much less system- 
size dependent even in the presence of a finite D c > 0. First 
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Figure 6: 2D cr(u>) for underdoped = 2/20, at different T/t and 
t'/t. 



the results for planar a(u>) in the underdoped case are pre- 
sented in Fig. 6. We note that the behavior is generally very 
similar to the one at intermediate doping in Fig. 2. Still, there 
appears to be already some build-up of a shoulder at, e.g., 
uj ~ 0.5 t for t' = 0, as a precursor of the mid-IR scale, hence 
a(uj) cannot be well described by a power law a oc uj~ v any 
more. Results for c^, = 0.1 can be considered as a crossover 
to substantially different behavior in the low-doping regime 
which we discuss furtheron. 

In Fig. 7 we present results for a(uj), as obtained for the 
largest system with N = 26 sites and a single hole Nh = 1. 
Very similar results were obtained also for Nh — 1 in systems 
with TV = 18, 20 sites. The behavior of <t(oj) is clearly differ- 
ent from the one belonging to the optimum doping in Fig. 2. 
The main difference comes from the emergence of a shoul- 
der at uj* which further depends on t! . While uj* ~ 0.5 t at 
t' = 0, it increases to uj* ~ t for t' = — 0.3i. The shoul- 
der starts to build up in <t(oj) for T < J which is a quite 
clear indication that the observed effect is related to the onset 
of short-range AFM correlations, confirming the relevance of 
the string-picture of the incoherent hole motion. 18 Lowering 
T < J we are in a regime of nearly constant a{uj < uj*) 
which at the same time induces the flattening or the saturation 
of p{T). It should be noted that our shoulder can plausibly 
be related to the mid-IR resonance in cupratesp*i£iA& in fact 
the position of both is even quantitatively well in agreement. 
At the same time a dip in u(uj) has also been observed in re- 
cent experiments in LSCO at low doping x = 0.08 and above 
the crossover temperature T > T* ~ 500 KfA£ as well as in 
YBCO^ 
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Figure 7: 2D a(ui) for low doping, calculated for Ch = 1/26 at 
different T/t and t'/t. 



Let us continue with the results for the normalized 2D re- 
sistivity ChP presented in Fig. 9 in a large range of tempera- 
tures T < t and for a wider regime of doping up to the 'op- 
timum' doping Ch = 0.15. For large T > 0.5 t all curves 
merge on the 'trivial' linear dependence dp/dT = (poks/cht 
with £ ~ 0.4, as follows, e.g., from the high-T expansion for 
p(T)p!2iil For the intermediate doping there appears a steady 
crossover at T ~ 0.25 t to the low-T linear law, as discussed 
in Sec. III. On the other hand, low-doping results reveal quite 
evidently a pseudo-saturation of p(T) within the T window 
T* < T < 0.5 t, closely related to the appearance of a flat 
and T-independent cr(u>) in this regime. We can compare our 
2D saturation value p sat with the estimate in RefsJ^iii, ob- 
tained by assuming a(ui < W) ~ const, with W being a 
typical band-width. Since (Hkin)/N ~ — 3.4te^(l — Ch) for 
2D t-J model 13 one gets from Eq. at low doping Ch — > 0, 



c^Psat ~ 0.37— po ~ 0.75p , (12) 



assuming an effective width of <j(o;) spectra W/i ~ 2 as esti- 
mated from Fig. 7, and po = h/ e o denotes the universal sheet 
resistivity. Our saturation value in Fig. 9 is well in agreement 
with the one following from Eq. il2\ . 



At decreasing T < T* < J our results indicate the on- 
set of a nondissipative contribution D C (T) > 0. Calcu- 
lated D C (T) are surprisingly consistent for different systems 
N = 18, 20, 26, as presented in Fig. 8. Plausibly one expects 
that in a large enough system additional scattering channels 
would lead to the broadening of the w = peak. Nevertheless 
its weight should not change considerably and also its width 
should remain at least narrower than the pseudogap scale oj*. 
Quite abrupt onset of D c > at T* — 0. 15 1 — 600 K is quite 
consistent with the experimentally observed scale T* at low 
doping, identified as the kink in p(T)*& To evaluate proper 
p(T < T*) we are missing the scattering mechanism within 
our small systems. For convenience we evaluate p(T < T*) 
by using spectra broadening e as elsewhere in our analysis. 
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Figure 8: Normalized charge stiffness D c / Ch as calculated for 
1 and N = 18, 20, 26 in the case of t' = 0. 




Figure 9: Normalized 2D resistivity Chp vs. T/t for different ch and 
t' = 0. 



Finally, we compare in Fig. 10 our results for Chp{T) at t' = 
with the experimental values for LSCO at low doping. 3 - 8 
First we note that the calculated values for Chp(T) in the low 
doping regime are larger compared to those at intermediate 
doping in Fig. 5 as observed in experiments. In the same way, 
it follows from Fig. 10 that Ch.p(T) is larger for Ch — 0.05 than 
for Ch = 0.1. Comparing our results for the lowest doping 
Ch = 1/20, 1/26 we can observe the behavior with a kink 
at T ~ T* similar to experiments. The values themselves 
are somewhat lower than those obtained by Ando et al£, and 
within a factor of 2 to previous data. 3 
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Figure 10: Normalized resistivity Chp(T) at low doping for t' = 
compared to experimental results for LSCO: x = 0.04, 0.07, crosses 
and squares, 3 and x — 0.04, 0.08, circles and triangles, 8 respectively. 



V. CONCLUSIONS 

Let us summarize and discuss our results for the optical 
conductivity u(ui) and the d.c. resistivity p(T) within the ex- 
tended t-J model. The TBC averaging method, used in our 
studies, provides a clear improvement of the FTLM exact- 
diagonalization studies of small systems. In particular, the 
advantage is evident in the low-doping regime, which is one 
of the central points in our analysis. 

Our study confirms that the transport properties substan- 
tially differ between the low-doping regime Ch < 0.1 and the 
intermediate doping, in our case c/, ~ 0.15. With respect 
to the NNN hopping t', our results indicate that the trans- 
port quantities are not much sensitive to its value, at least for 
T > Tf s which we are able to study numerically. General 
trend is that t' < 0, as appropriate for hole-doped cuprates, 
slightly increases p(T) as well as the shoulder (mid-IR) scale 
at low doping. At the same time, the characteristic behavior 
at intermediate doping moves with t' < towards the prop- 
erties at lower doping. In particular, properties at t' = and 
t' /t = —0.15 are qualitatively similar, while t' /t = —0.3 
already exhibits more pronounced deviation even at interme- 
diate doping. 

In this paper we present results for J/t — 0.3 only, as be- 
lieved to be the best parameter choice for cuprates. Still, addi- 
tional results obtained for J/t = 0.4 (not shown here) reveal 
that both o(uj) and p(T) do not change appreciably as far as 
we are dealing with the strong correlation regime J < t. 

At intermediate doping <j(uj) is well consistent with the 
anomalous scaling described within the MFL scenario. The 
overall behavior is close to the universal form, Eq. @, with 
C(uj) = Co in a wide range uj < u> c ~ 2 t and T < J. On the 
other hand, we can establish also deviations from the above 
simple form, in particular C(w) is an increasing function of 
uj. At uj > T the scaling is thus consistent with the power law 
a oc u>~ v with v < 1, in particular v ~ 0.65 for t! = very 
close to recent experimental result for BSCCO, 2 although one 
would expect better agreement with t 1 < results in the latter 
cuprates. 

The anomalous behavior of cf{uj) and C(uj) can be partly 
understood in terms of a simple representation of electron 



spectral functions, Eq. il It . where the quasiparticle damping 
T = y(\u\ + £T) is of the MFL form with large 7 > 0.3. 
It should be also noted that C(u>), Eq. il 11 . leads in general 
to an analytical dependence at small uj, i.e., Ta(uj — > 0) ~ 
B + A(u)/T) 2 rl in contrast to the simple C(uS) = C as- 
sumption. Still, we find that the same input for 7 which can 
describe the high-w behavior cannot quantitatively account for 
experimentally observed uj/T — > 0£ 

The origin of the anomalous MFL-like behavior of <j(u>) 
has not been clarified yet. The scaling range u> < 2 t is sur- 
prisingly broad, since one would expect that it might be also 
determined by J. On the other hand, the behavior qualita- 
tively changes for T > J. 1012 This is another sign, that we 
are dealing with an unusual frustrated fermionic system and 
possibly not with the quantum phase transition 2 which should 
possess only a single scale in the (uj, T) diagram. 

At low doping Ch < 0.1 the behavior changes qualitatively. 
cr(uj) exhibits a shoulder at uj* and flat region or a weak dip 
for uj < uj*, as found also in recent measurements in LSCO 15 
and YBCO. The shoulder is clearly related to the mid-IR 
peak found in experiments, whereby we find uj* > 2 J being 
in a reasonable quantitative agreement. The shoulder sets in 
for T < J and is plausibly related to the onset of short-range 
AFM correlations hence this confirms the magnetic origin of 
the mid-IR resonance and its spin-string interpretation>IiA& 
Another manifestation of the same phenomenon is the satu- 
ration of p(T) in the temperature window T* < T < J which 
is in agreement with experiments in cuprates. 14 

For T < T* we observe the onset of a coherent transport in 
our small systems, which shows up in a rather abrupt increase 
of the charge stiffness D C (T). The appearance of D C (T) > 
means that the scattering length L S (T) becomes larger than 
the system size y/N. In this connection, it is quite puzzling to 
explain why L s (T) is larger in the low doping regime than at 
intermediate doping. I.e., we find D c ~ down to T ~ T/ s 
at = 0.15, whereas evidently we get D c > at the same 
T in the low-doping regime. This apparent contradiction is, 
however, consistent with recent experimental confirmation of 
a quite coherent transport at low doping for T < T* and the 
concept of a nodal metal£A& 

While qualitative behavior of p(T) is well in agreement 
with experiments, let us finally comment on the quantitative 
comparison. At intermediate doping the normalized resistivity 
Ch,p(T) lies tolerably within the range of existing experimen- 
tal values for LSCO, taking into account also some variation 
between different published results. At low doping our results 
seem to underestimate somewhat p(T). Still, the discrepancy 
(depending again on different experimental data) is at most of 
a factor of 2. Whether this rather modest disagreement is an 
indication of the possible relevant role of degrees of freedom 
outside the t-J model, remains an open question. 
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